From Rotating Atomic Rings to Quantum Hall States 
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Considerable efforts are currently devoted to the preparation of ultracold neutral atoms in the emblematic 
strongly correlated quantum Hall regime. The routes followed so far essentially rely on thermodynamics, i.e. 
imposing the proper Hamiltonian and cooling the system towards its ground state. In rapidly rotating 2D har- 
monic traps the role of the transverse magnetic field is played by the angular velocity. For particle numbers 
significantly larger than unity, the required angular momentum is very large and it can be obtained only for spin- 
ning frequencies extremely near to the deconfinement limit; consequently, the required control on experimental 
parameters turns out to be far too stringent. Here we propose to follow instead a dynamic path starting from 
the gas confined in a rotating ring. The large moment of inertia of the fluid facilitates the access to states with 
a large angular momentum, corresponding to a giant vortex. The initial ring-shaped trapping potential is then 
adiabatically transformed into a harmonic confinement, which brings the interacting atomic gas in the desired 
quantum Hall regime. We provide clear numerical evidence that for a relatively broad range of initial angular 
frequencies, the giant vortex state is adiabatically connected to the bosonic v = 1/2 Laughlin state, and we 
discuss the scaling to many particles. 



PACS numbers: 73.43.-f, 05.30.Jp, 03.75.Kk 

While coherence between atoms finds its realization in 
Bose-Einstein condensates [ lJO, quantum Hall states (H are 
emblematic representatives of the strongly correlated regime. 
The fractional quantum Hall effect (FQHE) has been discov- 
ered in the early 1980s by applying a transverse magnetic 
field to a two-dimensional (2D) electron gas confined in semi- 
conductor heteroj unctions 0. Since then, FQHE has never 
stopped to intrigue the scientific community due to non-trivial 
transport properties and exotic topological quantum phases 
l6l . Such interest has also influenced the research in ultra- 
cold atomic gases, which in the last decade have been success- 
fully exploited as a highly controllable playground for quan- 
tum simulations of many -body physics 0. The large versa- 
tility of these setups allows one to confine atoms in 2D har- 
monic traps and to impose an effective magnetic field either 
by rapid rotation 00 or by laser-induced geometric gauge 
potentials [ 9 1. In principle, such opportunity should allow one 
to experimentally explore the bosonic version of QHE, even if 
unfortunately it has been hitherto elusive. 

From a theoretical point of view, a variety of interesting 
ground states (GS's) have been identified for Bose gases as a 
function of the effective magnetic field (7] [8) . At zero field, 
i.e., without rotation, the particles undergo Bose-Einstein 
condensation [1, 2 ] and the atomic ensemble is superfluid. 
Differently from a rigid body, a superfluid of N particles reacts 
to rotation with the formation of quantized vortices, whose 
number N$ increases with the rotation frequency. At large fill- 
ing factor v=A/yAfy,>10an ordered vortex lattice is formed. 
For v < 10, the lattice melts because of quantum fluctuations, 
which signals the breakdown of the mean-field description 
and the access into the FQHE regime. The filling factor is 
now better defined as v = N/m mSLX with m max the maximum 
angular momentum occupied by single particles. FQHE states 
are obtained for values of v of order unity, which correspond 
to very large total angular momenta LocN 2 . 



Like in solid-state physics, most of the preparation pro- 
cedures employed so far in rotating atomic ensembles ap- 
proached the GS by cooling down the system with a fixed 
Hamiltonian. By contrast we explore in this paper an alter- 
native method that consists in starting from an easily prepara- 
ble state (typically uncorrected), following a dynamic route 
by changing an external parameter, and eventually obtaining 
the desired state. This strategy has been successful for the 
experimental investigation of the superfluid to Mott insulator 
transition in optical lattices 1 10 ]. We propose to implement it 
to reach Quantum Hall states with the following steps: (i) We 
engineer a Mexican-hat trapping potential by superposing a 
standard harmonic trap with the repulsive potential created by 
a "plug" laser beam, which is focused at the center of the trap, 
(ii) By stirring the gas, we prepare the N bosonic atoms in 
a giant vortex state, corresponding to the lowest energy state 
of the Mexican-hat potential for a given angular momentum 
L. (iii) The stirring is removed and the plug is adiabatically 
switched off. (iv) In the final harmonic trap, we obtain the GS 
with the initially imparted angular momentum L, thanks to ro- 
tational symmetry. We show that if L = N(N — 1) then the 
2-body contact interactions drive the gas into the celebrated 
bosonic V = 1/2 Laughlin state HI . 

Single particle physics 

In experiments with rotating atomic gases, particles are 
usually trapped by a harmonic potential and stirred by time- 
varying magnetic field or auxiliary laser beams Q. In the 
frame rotating at angular speed H = £2z, the Hamiltonian of a 
single particle in the harmonic trap of frequencies (ft), (0, co z ) 
can be written as 

^a P = ^ (P - A) 2 + y {co 2 - tf){x 2 +y 2 ) + ™co 2 z z 2 , (1) 
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Figure 1: Mexican hat potential:(a) Bosonic atoms are confined in a combined trap with (i) an isotropic harmonic confinement and (ii) the 
dipole potential created by a blue-detuned, gaussian laser beam that plugs the trap center and pushes the particles away from this point. The 
resulting potential exhibits a Mexican-hat shape, (b) Under fast rotation, the single-particle energy spectrum exhibits a Landau Level picture 
(red dashes), where levels are arranged in quasi-degenerate manifolds as in the case of a purely harmonic trapping (blue asterisks). The level 
plot has been drawn for 5 = 0.09, a = 3.0 and w 2 = 8.0, parameters that we will use in the many-body problem for N = 9. The single-particle 
minimum of energy can be adjusted to any desired value on angular momentum by tuning the waist and the power of the plug beam (m = 8 in 
the present case), (c) In order to enter FQHE regime of the many-body interacting system, we propose to switch off the plug beam, eventually 
recovering the usual parabolic form. 



with A = MCI x r = M£l(—y,x,0). In the following we sup- 
pose that all relevant energies are much smaller than hco z , so 
that the motion along the z direction is frozen and the prob- 
lem is effectively two-dimensional. In the limit of centrifu- 
gal deconfinement £1 — » CO, the system is formally equiva- 
lent to bosons of charge q = 1 in uniform magnetic field 
B = V x A = 2M€lz. From now on, we express energies and 
lengths in units of kco and ^k/Mco, respectively. It is well 
known [7 ] that the problem can be rewritten as two decou- 
pled harmonic oscillators ^tmp = (2cft a + 1) + S(b^b — eft a), 
in terms of ladder operators a, b, and 5 = 1 — £l/co is the fre- 
quency offset. Every state is labeled by the occupation num- 
ber n a , rib of the two modes, and it is denoted as | Yn b ,n a ) • Note 
that a gauge field similar to the one entering into ([T]) can also 
be induced by geometric phases instead of rotation [ 9 ]. The 
scheme outlined in the present paper should work equally well 
in this case, the only significant difference being that ([T]) is 
now the single-particle Hamiltonian in the laboratory frame, 
instead of the rotating frame. 

In the limit <5 <C 1, the quantum number n a identifies dif- 
ferent manifolds called Landau Levels (LL). Within each LL, 
the states (labeled by rib) are quasi-degenerate due to the small 
separation energy 5. The quantity m — rib — n a is the angular 
momentum of the particle. In the lowest Landau level (LLL), 
n a = and the one-body eigenf unctions assume the simple 
expression \l/ m o(z) = -i=z m e~^ / 2 , where z now denotes 
the position in the complex plane (z = x + iy) , with energies 
E m = mS and angular momentum m. 

The first key feature of our proposition is to replace the 
ordinary harmonic potential with a Mexican-hat one, like in 



Fig[T] This can be done by shining the center of the harmonic 
trap with a laser beam prepared in a circular, Gaussian TEMqo 
mode 1 1 1 ]. When the laser frequency is chosen larger than the 
atomic resonance frequency ('blue detuning'), the laser beam 
creates a repulsive dipole potential proportional to the light 
intensity. The beam is chosen to be perpendicular to the xy 
plane and the dipole potential is of the form 

U w (x,y) = aexp [-2(x 2 +/)/w 2 ] , (2) 

where w is the laser waist and a is proportional to the laser 
intensity. The sum of the harmonic potential (x 2 +y 2 )/2 and 
of U w (x,y) has a bump in x = y = in the laboratory frame 
when a > w 2 /4. 

At moderate intensities of the plug, as the ones employed 
in our preparation scheme, the classification of single-particle 
energy eigenstates in terms of LL remains valid (see Fig. [T] 
and Methods). In the LLL the single-body energies are in 
good approximation: 

/ 2 \~( m+1 ) 
e m = mS + a f 1 + ~2 J • O) 

At fixed laser parameters a and w, the angular momentum 
m that minimizes £ m is a decreasing function of the rotation 
frequency offset 8. We denote by S m the value for which the 
level crossing £ m+ i = £ m occurs. The LLL state with angular 
momentum m is thus the lowest energy state when 5 is chosen 
in the interval S m < 8 < <5 m _i, whose width is 

, 2 x2 , 2 x -(m+2) 
hn = S m -\ - S m = a ( — 2 I + — 2 J 
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Later on, we will be interested in choosing a specific value 
m = t and in maximizing the width of the stability window. 
This can be done, at fixed intensity a, by choosing w 2 = t. 
The central rotation frequency in the stability window for £ 
then corresponds to 



1 



(S £ _ 1 + S £ ) = 2a 
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(4) 



Notice that for large values of we get 8^ ^ ai 1 , thus if we 
want to keep it sizable, we have to choose a^c £. 



Many-body physics 

In the context of cold bosonic gases in the LLL subspace, 
two-particles interactions can be modelled by the contact po- 
tential 



(5) 



whose strength is given by the adimensional parameter C2 = 
\^7la s /a z , where a s is the 3D s-wave scattering length and 
a z = ^h/Mco z is the size of the ground state in the strongly 
confined direction [7 , 8 ]. Within the kernel of J%2 the V = 1/2 
Laughlin state 



^Lau = 



(6) 



has the lowest total angular momentum = N{N — 1), or 
equivalently the angular momentum per particle ^Lau =N—1. 

We first recall the practical difficulties to attain the FQHE 
regime via a thermodynamic route for a pure harmonic con- 
finement in the xy plane. In a typical experiment with 87 Rb 
atoms (a s — 5nm), a longitudinal frequency (O z /2n = 50kHz 
gives a z = 50 nm and an interaction parameter C2 = 0.5. Then, 
already for a modest number of particles, the Laughlin state 
is reached only for rotation frequencies £1 extremely close to 
the centrifugal limit CO. For N — 9 and C2 = 0.5, we find us- 
ing exact numerical diagonalization that the Laughlin state is 
the GS only for <5L au < 5.5 10 -3 (see figure]^). This very 
low threshold makes it difficult to transfer the desired angular 
momentum to the gas. Indeed when the stirrer consists in a 
rotating anisotropic potential e(x 2 — y 2 )/2, the corresponding 
anisotropy £ must be chosen smaller than 5 to avoid a dynami- 
cal instability of the center-of-mass motion [12]. One has thus 
to restrict to extremely weak stirrers, with £ in the 10 -3 range. 
However one must also choose £^> u, where u is the strength 
of the static anisotropy defect u(x 2 — y 2 )/2, otherwise the gas 
cannot be effectively set in rotation. Unfortunately, in real- 
istic traps the typical values of u are at least of the order of 
10 -3 . Consequently it is quite problematic to fulfil simulta- 
neously these various constraints. Lastly, we mention that in 
such a thermodynamical procedure the temperature has to be 
kept below 5, i.e the gap between the Laughlin state and the 



GS with L = LLau + 1 (the interaction energy is zero in both 
cases), which is a very stringent requirement. 

Here we propose a different point of view where the FQHE 
regime can be tackled from a dynamical perspective, with a 
two-step procedure. The first step is addressed in this section 
and it consists in the preparation of a giant vortex state of N 
particles in the Mexican-hat potential of Fig. [T^, with the de- 
sired angular momentum L = L^ au . This preparation is easier 
than the direct production of the Laughlin state, thanks to the 
favourable parameter sensitivity of the Mexican-hat potential. 
More precisely the increased moment of inertia of the gas en- 
ables one to reach L oc N 2 in a relatively broad interval of £1. 
The second step involves the adiabatic transformation of the 
giant vortex state in the Laughlin state, and it will be analyzed 
in the next section. 

In the limit case of no interactions, every single particle in 
the Mexican hat potential should rotate at angular momentum 
t — L/N. The bosons then condense in the giant vortex state 



N 
i=l 



(7) 



similar to those already observed in f\3l and theoretically an- 
alyzed in [14]. As shown in the previous section, the window 
of stability /[^Lau] is optimized for w 2 = ^Lau- For N = 9 the 
constraint a > w 2 /4 imposes a > 2. We choose in the fol- 
lowing a = 3, which leads to 5L au e (0.081, 0.101). 

In the presence of interactions, the interval of values for 5 
leading to a ground state with L — can be determined 
either from a Bogoliubov analysis or from exact diagonal- 
ization. The main role of the interactions is to deplete the 
contribution of the mean angular momentum i in favour of 
neighbouring ones £±q, with q <C L The Bogoliubov anal- 
ysis is well suited for strong plugs (a > 1) where the de- 
pletion is small, while for intermediate regimes a full many- 
body numerical treatment is needed (see Methods). Thanks 
to the angular momentum conservation, the exact diagonal- 
ization can be performed in each L sector separately and the 
conjugate variable 5 simply yields the energy shifts L<5. The 
phase diagram as a function of (a, 5) is presented in Fig|2j) 
for N = 9, C2 = 0.5. It strongly supports the sketch drawn 
before for non-interacting particles. In particular the com- 
puted ground state for a Mexican-hat potential with a = 3 
possesses the required angular momentum L^ au for the inter- 
val <5L au G (0.084, 0.105), very close to the one in absence of 
interaction. This corresponds to a ~ 10% frequency difference 
between £1 and 0), which is notably larger than the typical stir- 
rer anisotropy needed to set a gas in rotation. This ensures that 
the preparation of the giant vortex state with L = L^ m should 
be rather robust. 

It is important to stress that the Mexican-hat potential is 
employed just for the scope of injecting the right quantity of 
angular momentum L^ m , and not for producing the Laughlin 
state itself. The situation considered here is thus completely 
different from former proposals suggesting to find a tradeoff 
between £1 and a that optimizes the fidelity with the Laughlin 
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Figure 2: Phase diagram for N = 9 and C2 = 0.5. Exact diagonal- 
ization in a truncated LLL basis (m <2N = 18) is performed sepa- 
rately for each sector of L, thanks to rotational invariance. Energies 
are then shifted by the total angular momentum term Ld to draw 
the phase boundaries, (a) In absence of the plug beam, a = 0, the 
window of stability for Laughlin angular momentum (LLau — 72) is 
narrow and extremely close to deconfinement limit, (b) Conversely, 
the region with L]^ m opens up and drifts away from 8 = as the plug 
power a is ramped up (at constant w 2 = 8); the same happens for 
other large angular momenta around it. For large values of a (typi- 
cally larger than unity), the GS found when varying 5 are essentially 
non-correlated states, where all atoms accumulate in the same giant 
vortex state. Consequently the total angular momentum L undergoes 
jumps of size N, corresponding to the addition of one flux quantum 
to each particle. 

state fT5l . 



Adiabatic evolution 

Once the gas has gained the desired angular momentum 
L = N(N — 1) via equilibrating in the giant vortex state |7]), 
the stirrer at frequency £1 can be suppressed. The situation 
becomes rotationally symmetric and the total angular momen- 
tum is thus conserved. Then, the slow removal of the laser 
plug will result in a redistribution of particles around the mean 
angular momentum by repulsive interactions. Such a re- 
distribution reaches a paradigmatic form in the unplugged har- 
monic trap, where the Laughlin state ^ has no interaction en- 
ergy anymore. From a technical point of view, we note that in 
absence of stirring we can now look for the GS of the gas in 
the laboratory frame, within the subspace of the L z that had 
been imparted to the cloud during the stirring phase. 

The system will follow the instantaneous GS if the un- 
plugging path can be followed slowly enough to satisfy the 
adiabatic condition | (¥ | (dJt?/dt) | <C Aj, where j 
represents an excited eigenstate of energy Ej of the instan- 
taneous Hamiltonian Jt? 9 and where Aj = Ej — Eo fT6l . We 
have checked numerically that the most stringent constraint 
originates from the first excited state *Fi and we thus focus 
our discussion on this state. Once the gap A = E\ — Eq is 
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Figure 3: Finite size scaling analysis. Data collapse of (a) the 
gap and (b) interaction energy as a function of a for an ini- 

tial value of ji = C2N 1 / 2 / (2k 3 / 2 ) ~ 0.15, corresponding to C2 = 0.5 
for N = 9. The inset in plot (a) shows the finite size scaling of the 
Laughlin gap in the harmonic case (a = 0). Our estimates gives 
^Lau — 0.097(1) in the thermodynamical limit, (c) Scaling of the 
function F a = A -2 \( x ¥q\ (dU w /da) whose integral provides 
an estimation for the adiabatic time (inset). At sizable enough a's 
(> 0. IN) the curves collapse from above, giving a total time T which 
is approximately constant with N ("x" points). Alas the presence of 
a pronounced bump for small a's leads to T °< Af ("+" points). Al- 
ternative strategies that lead to lower adiabatic times are discussed in 
Fig- 

known from exact diagonalization, the changing rate of laser 
intensity a and/or its rescaled cross section a = w 2 /(N — 1) 
with time t is determined by the condition 

F a da + F G do<z:dt , (8) 

where F x = A -2 | (^o I (dU w /dx) |¥i)| is the matrix element 
of the potential variation in x = a, a. The minimal total time 
T required for adiabaticity will then be the integral of those 
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functions along the chosen path, T = J (F a da + F G do). In 
the following, we first consider the case where the waist is 
kept fixed, which is experimentally straightforward since it 
involves only a variation of the laser intensity; then we address 
the general case of changing of both a and a. 

We have performed numerical simulations for up to N = 10 
particles, in a LLL truncated single-particle basis m < 2N, in 
order to test the validity of the adiabatic approximation (see 
Methods). For the chosen test case of N = 9, C2 = 0.5, ramp- 
ing down the intensity from the initial value a = 4.5 at con- 
stant a = 1, we obtain T « 43 (in units of co~ l ). Such a value 
of T is a reasonable time in state-of-art experiments, estab- 
lishing the feasibility of our scheme for N = 9, as opposed to 
the procedure involving a purely harmonic rotating trap. 

The exponentially increasing dimension of the Hilbert 
space and the strong correlations involved ward off going 
much further than N = 10 particles with the exact diagonal- 
ization method. To infer the behavior of larger samples, we 
performed finite size scaling of the relevant energies using 
the Bogoliubov approximation (see Methods). Our scheme 
requires the preparation of the gas in the ring with a oc N 
and a — 1, for which the chemical potential goes as ji ~ 
C2N 1 / 2 / (2n 3 / 2 ). The LLL approximation requires jl < 2, and 
working at fixed /i implies C2 N~ 1 / 2 . We then deduce that 
the energy gap behaves like A oc N~ 1 / 2 and the interaction en- 
ergy as (J#2) °* N (see Methods). We have plotted in FigJ^a)- 
(b) the variations of the gap A and the interaction energy (J#2) 
using C2 = 1.57V -1 / 2 . The expected data collapse is well veri- 
fied for values of a/N larger than 0.1. 

A finite-size scaling can be performed also for the quantity 
F a entering the adiabatic condition ^ and we plot the result 
at fixed waist a in Fig|3j:. This function takes its largest values 
in the interval a G (0, a c ) with a c — 0.1 /A/. We can therefore 
decompose the adiabatic path into two successive parts. In 
the first part the plug laser intensity a is decreased from Of = 
0.5 /N down to a c and this can be down adiabatically in a 
relatively short time T\ = J^F a da ~ 5, independent of the 
number of particles (inset of Fig |5j:). In the second part (0 < 
a < a c ) the breakdown of the scaling A oc A/ -1 / 2 imposes a 
slowdown in the reduction of the plug intensity. The time T2 = 
Jq c F a da needed for this second part actually show a linear 
increase with N, hampering the feasibility for more than a few 
tens of bosons. 

The sequence of spatial density profiles during the time 
evolution is depicted in Fig|4^. From such a sequence it is 
clear that the gas starts feeding the trap center at the end of 
the path. An alternative strategy to ramp down the plug con- 
sists on reducing its waist w while maintaining a constant, 
this being performed in practice (up to the diffraction limit) 
using a motorized focusing optical element. The correspond- 
ing evolution of the density profile in the trap is represented in 
Fig|4j3 and in our specific case it can be covered adiabatically 
in half time with respect to the above situation (see Fig|5]). 

A natural extension of our analysis is to consider a simulta- 
neous ramping of a and <7, in order to minimize the total evo- 
lution time while fulfilling the adiabaticity criterion. To this 



aim, constrained optimization techniques can be implemented 
using the data of the vector (F G ,F a ), represented in Fig [5] 
Experimentally, another effective way of reducing the adia- 
batic ramp time is to increase the interaction coupling con- 
stant C2, hence the gap, via either Feshbach resonances fTTl or 
a tighter longitudinal confinement co z . For a ramp of a only, 
our numerical calculations with N — 9 give T ~ 65,43, 20 for 
C2 = 0.33,0.5, 1.0, respectively, corresponding to the empiri- 
cal scaling law T « 20 . 

Finally we briefly address the consequences of some of 
the unavoidable experimental imperfections on the proposed 
scheme. The two principal perturbations that we can fore- 
see are the imperfect centering of the plug beam and the 
residual trap anisotropy. We model these defects by writ- 
ing the dipole potential created by the plug beam as U w = 
aexp [— 2[(x — v) 2 + y 2 ]/w 2 ] , and by adding the term u(x 2 — 
y 2 )/2 to the single-particle Hamiltonian to account for the 
static anistropic defect. Here v and u are dimensionless coeffi- 
cients characterising these imperfections. These two coupling 
terms break the rotation symmetry: in their presence, the an- 
gular momentum is not a conserved quantity anymore and the 
gas will undergo a cascade from L = down to states with 
no angular momentum, by populating the first excited LL. To 
get a conservative estimate, we impose the very stringent con- 
dition that the total angular momentum remains unchanged 
over the adiabatic ramp time, and we estimate the correspond- 
ing constraint on u and v using time-dependent perturbation 
theory (see Methods). The constraint on u is certainly the 
most challenging one. We find that the maximal tolerable trap 
anisotropy u max £ 2AL au /Af ~ O.2C2/N. Taking u ~ 10 -3 as 
a realistic trap anisotropy, we find that our scheme should be 
operational for atom numbers up to Af max = 100 for C2 = 0.5. 

Detection 

One of the simplest techniques to probe cold atomic setups 
consists in taking time-of-flight (TOF) pictures |3|. The ab- 
sorption image of the density profile expanded after releas- 
ing the harmonic confinement contains indeed useful informa- 
tions about the initial situation in the trap. In the specific case 
of bosons in the LLL regime, the density profile is self-similar 
in time and the TOF picture simply magnifies the original par- 
ticle distribution in the trap fT8l . Given the direct connection 
between single-particle angular momenta and orbital radius 
(see Methods), a TOF image allows one to compute the angu- 
lar momentum. The V = 1 /2 Laughlin state with N particles 
exhibits a fairly flat profile of density 0.5 inside a rim of radius 
~ y/N. Observing such TOF images would be already a first 
hint that one has effectively reached the QHE regime. 

Multi-particle correlations offer even more insight into the 
many-body state. These correlations are directly accessible 
if one uses a detection scheme that can resolve individual 
atoms with sub-micron resolution ||T9j [20). Alternatively the 
two-body correlation function can be tested at short distances 
using the resonant photo-association of spatially close pairs 
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Figure 4: Density profile during adiabatic evolution (N = 9). The leftmost panel corresponds to a giant vortex like structure, whereas 
the rightmost one depicts the flat disk shaped profile of the Laughlin state. In the upper row cr = 1 is kept constant while a = 
1., 0.5, 0.4, 0.3, 0.2, 0.1, 0. The last part of the ramp down procedure < a < 0.1 is the slowest, due to the large value of F a in this re- 
gion (see Fig. [5J:). In the lower row we squeeze the laser waist a = 1., 0.5, 0.25, 0.125, 0.025, 0.00625, 0. at fixed intensity a = 1.: particles 
spread towards the inner part of the trap in a different way, corresponding in a lower value of F a and faster allowed rates of change. For 
systems within LLL, density profiles after trap release and time-of-flight imaging will simply display rescalings of these pictures. 
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of initial 5. 

Production, observation and control of anyons is one of the 
most intriguing quests motivating the considerable efforts to- 
wards FQHE regime. Anyons are quasiparticles with the pe- 
culiar property of satisfying non common braiding rules when 
moving around each other. In the Laughlin case the anyonic 
excitations are quasiholes Y\i(zi ~ -ZqiO^Lau* which can be pro- 
duced and controlled by impinching a narrow strongly repul- 
sive laser beam at position z q h as put forward in [ 22 ]. A further 
feature of our proposal is that addressing a giant vortex with 
£ > N permits in principle to get a final state with a whole 
manifold of anyonic quasiholes and to study its exotic proper- 
ties. 



Figure 5: Map of adiabaticity requirements. Absolute value of the 
vector (F(j, F a ) is plotted in the coloured map for N = 9, evidenc- 
ing the large value of F a at large a = w 2 /(N — 1) and small a, as 
well as the more favorable condition if one uses a reduction in time 
of the beam waist. The two paths described in the text give T ~ 40 
(solid blue line) and T ~ 20 (dashed blue line). Superimposed white 
arrows represent the directions of the vector (F a ,F a ). This plot can 
serve for conceiving more intricate paths with the help of optimiza- 
tion techniques. 



I2TII . The amount of produced molecules is indeed directly re- 
lated to the correlation function g^ (0), which is also in direct 
correspondence with the interaction energy (J^)/c2 (Fig|5J)). 
Since the Laughlin state belongs to the kernel of Jf%, its pres- 
ence will be signaled by a strong suppression of two-body 
losses. Moreover, in a strict analogy with solid state physics, 
we can imagine an experiment to measure FQHE plateaus in 
physical quantities. Namely, by varying the rotational offset 5 
in the giant vortex preparation stage it is possible to change L 
by steps of N, i.e. move the penetrating magnetic flux in units 
of single quanta. Removing now the plug, the system will fall 
in a sequence of incompressible FQHE states: the final (0) 
is expected to display plateaus at discrete values as a function 
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Methods 

LLL approximation. Freezing the longitudinal degrees of free- 
dom by a large (D z , the single -particle Hamiltonian ([TJ can be rewrit- 
ten in dimensionless units as 

— a* a), 

with = —d z + z/2 and = — d z +z/2 in terms of z = (x+iy). 
Eigenstates with angular momenta I = m — p and energies £ m ^ p = 
2p + Id are built by iteratively applying the ladder operators 



7 



and on the vacuum i//q,o = (z|0, 0) = e \ z \ 2 / 2 /y/n, i.e. \m,p) 



(a )'_(P) |q q^ -phe explicit wavefunctions are 



jp\m\ 



m,p ■ 



Vo,o- 



-5 y (- 

to 



qj \q 



qi\z 



|2(P-«) J 



which can be rewritten in terms of confluent hypergeometric func- 
tions U(—p,m — p+ 1, |z| 2 ). Energy levels for small 5's are or- 
ganized in quasi-degenerate manifolds called Landau Levels (LL), 
labelled by the integer p and separated by an energy gap 2. When 
dealing with many-body problems, the usual approximation is to cut 
down the single-particle Hilbert space to the lowest LL (LLL) p = 
where wavefunctions (apart from Gaussian weight) are analytical in 
z, being \p m § = x^q q • z m j\fm\. This is well justified and valid if 
the chemical potential remains well under the LL gap value 2; fur- 
thermore, the LL mixing due to the interaction term in the strongly 
correlated FQHE regime is negligible ['23 1. 

Here we show that the LLL approximation remains valid even in 
the presence of a plug laser (2}. The matrix elements within the LLL 
are 

2 \ ~(m+l) 

(9) 



(m,0\U w (r)\m,0) = a[ 1 + 



which provide the energy shifts due to U w at the first order of per- 
turbation theory. The rotationally symmetric U w only couples states 
with the same angular momentum and different LL labels: 

(m + p,p\U w \m + p,p) ~ (m,0\U w \m,0) - + ^ 

(m,0\U w \m + p,p) ~ {y/m/y^Y (m,0 m,0) . (10) 

Due to the exponential decay with m in Eq.{9}, the plug laser affects 
mainly low angular momenta, localized inside a circular area of ra- 
dius ^fm < w. For large w the energy shift of higher LL is quite 
similar to the LLL one, thus almost preserving the distance 2 be- 
tween adjacent LLs. The inter-LL terms in the illuminated region 
m < w 2 are reduced by the prefator ^/mw~ 2 , varying exponentially 
with the LL index p. Formally, one is allowed to use LLL approxi- 
mation only for those m such that the matrix elements {To} are <C 2. 
In practice among the exact eigenstates m? o of ^ rap + U w , those 
that differ significantly from the unperturbed states \j/ m ^ have a very 
small occupation in the many-body solution: the global LL mixing 
Lm[(l — l(V / m,o|0m,o)| 2 )(^m)] is bounded by 1% for all the simula- 
tions presented (see Fig[T]:). 

Numerics. LLL approximation reduces the particles' degrees of 
freedom to one, the angular momentum m, and the many body sys- 
tem is then described by the Fock basis \non\ ...). The Hilbert space 
available for a single particle is further cut to < m < 2N in or- 
der to realize numerics without affecting the correct description of 
the Laughlin state and its lowest excitations; occupation number in 
higher m's never exceeds negligible amounts. Even with this strong 
reduction, the dimension of the largest Hilbert subspace considered 
for N particles grows as - 1.75 • 10°- 74A ^ 2 (~ 4 • 10 5 for N = 10). 
Within the LLL approximation, central contact interactions of the 
form ([5} could be written in terms of a single Haldane pseudopoten- 
tial(24]: 



Here we introduced the second quantization operator dm, p , which 
creates a particle in V6n,p» an d used the simplified notation = 



for LLL states. The coefficients in (TT) account for pairs having total 
angular momentum mo with null component in the center of mass 
frame. This leads to a sparse matrix form for the interactions, with an 
average filling per row growing as ~ 0.2-N 292 . For N = 10 particles 
and L = Ll eu = 90, we need ~ 1 Gb RAM to store the Hamiltonians 
and less than one hour CPU-time on a single-core 3GHz desktop 
processor to diagonalize a single instance of the problem. 

Condensate in the ring. Given the noninteracting energies (3} 
with a minimum in m = £, the GS of M\ = )Z m £mdmd m is given by 
the giant vortex (7}. The angular momentum £ gets depleted by the 
insertion of two body interactions {TT} in favour of the nearest ones 
£ + q, £ — q, with q <C L For a condensate of Nq particles in m = £, 
the most dominant terms in J%2 are 



2^ 



C2 



N 2 +N £ 2p}p q + 2pl q P- q + p l P- q + pl q $ 



q^0 



where new operators P q = d^ +q have been defined and the Stirling 

approximation n\ « y/2jie~ n n n+ ^ employed. The coupling gets 
renormalized by £~ 1 / 2 as a consequence of the wavefunction local- 
ization on a ring of length 2%\T£- Eliminating No by the number 
operator N = Nq + PqPq, the overall Hamiltonian M{ + reads 



= Net + 8 -N 2 + ^(e q +gN)f$lfi q 



L 1 



(12) 



with g - 



: ££+ q — ££~ 2aq 2 / (e£) 2 for small q and e is 
the exponential constant. Under the proposed scaling a oc£ = N— 1, 
increasing N enhances the importance of interactions with respect to 
single-body energies. 

The quadratic bosonic model (T2| can be exactly solved |25 | by 
the Bogoliubov transformation P q = u q r) q — v q T\L q with u 2 — v 2 = l, 
by which it reads 



3? = Ne i + g -N 2 + L -Y,{\ 



1 



'^0 



-£ q -gN)+ £A^*V 

q^0 



Due to the positiveness of A q , the GS |3>o) (an approximation to 
the exact GS |¥q)) is given by the vacuum of quasiparticle ex- 
citations, f^l^o) — 0» V#. The quasiparticle spectrum is = 

N) 2 — g 2 N 2 and the occupation number of the state q is 

1^. The depletion is the fraction of particles out- 



l ( £ q +gN 



y 4 - 2 \ A q 
side the condensate 



A^o 1 v 2 

q^0 



1 

2N 



q^0 



q 2 + Q 2 



_\]{q 2 + Q 2 ?-Q A 



(13) 



with Q 2 = gN(eX) 2 /2a. The expression for v 2 converges as q~ 4 
for high momenta, and no infrared divergencies appear since angular 
momenta are quantized in integers. For large g, Eq.|T3|) becomes 
1—No/Ntt Q\ogQ/{y/lN) which vanishes foxN^°° only if Q oc 
N l ~^, with f > 0. In the same limit the GS energy turns out to be 

gN gQ g 
'2 3V2 4 



E^N[e £ 
and the chemical potential 

dE m 



e (gN 



1/2 



(14) 
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whose leading term is the expected result from the Gross-Pitaevskii 
approach in the LLL. Working at constant ji < 2, as required by LLL 
approximation, implies a scaling C2 = CAf -1 / 2 , with C = iVn^ji 
(see Figji]). Moreover, a N implies that Q oc A/ 4 / 2 , ensuring the 
vanishing of the depletion fraction. The radial confinement is rather 
strong, since the standard deviation Aq grows only as Af 1 / 4 , as de- 
duced from the calculation (Aq) 2 = Zqfvjj/N na V2Q 3 /(3N). 

The energy gap to the first excitation |4>i) = Tjjljlj |4>o) with 
the same L reads 

A = A 1+ A_ 1= 2giV^(e-2 + l)2-l«^ 

and vanishes as N~ 1 / 2 just as C2, i.e. the energy scale of the final 
Laughlin state (Fig[3j. Within such Bogoliubov analysis, it is also 
possible to determine the scaling of many other interesting quantities: 
e.g. the interaction energy (J%) = C2dE/dc2 scales as N. 

Robustness against trap defects. The main experimental de- 
fects that may hinder our protocol are a residual static quadrupole 
anisotropy u and an off-centering v of the plug beam, which are 
described by the single-particle potentials = u(x 2 — y 2 )/2 and 
J% = aexp [— 2[(x — v) 2 + y 2 ]/w 2 ] — U w , respectively. Both terms 
break the rotation symmetry and couple manifolds corresponding 
to different total angular momenta. We consider first the cou- 
pling Jif u since it turns out to have the largest impact for prac- 
tical conditions. Its second-quantized expression can be written 

JP« = M° ] + M l) + M 2) with 

^ = - £ y/(m+ l)(m + 2) (4 +2 / m ,, + H.c.) , 

m,p 

■k w = |LV^+i)(4-i,p+i^+ H - c -) . 

m,p 

^ = \ 1 V(P+ l)(P + 2) (< P+2 ^ )P + H.c.) , 

m,p 

where creates a particle in the one-body state \\f m , p . In the ab- 
sence of trap defects, we assume that the adiabatic approximation 
is valid and that the system is at any time in the instantaneous GS 
|*Fo) °f me time-dependent Hamiltonian ^g, given by the trapping 
potential plus the interaction energy. The state l^o) has angular mo- 
mentum L, energy Eq and belongs to the LLL (p = 0). We estimate 
the influence of Jf? u by calculating perturbatively the depletion of 
the probability for finding the system in l^o)- The states that can 
be reached are given at first order by the action of the various terms 

in J4? u on l^o)- The first contribution induces a coupling to 

other LLL states having angular momentum L± 2. The contributions 

and connect l^o) with states having both angular mo- 

mentum L — 2 and belonging to the first and second excited LL, re- 
spectively. We denote these normalized states as 00 l^o) 

and^V^Vo)- 

In analogy with the textbook problem of Rabi oscillations, we 
find that the initial state l^o) is protected from depletion towards 
another eigenstate \^a} of J% with energy E a if the energy detun- 
ing \E a -Eo \ is much bigger than the coupling T a = |(¥ a |^i|¥o)|- 
In the opposite case where T a ^> \E a —Eq\, the population of I^Pq) 
slims down as 1 — sin 2 (T a t). In such a resonant case, we should 
have r a « I 1 " 1 where T is the total time of evolution, to avoid any 
significant leakage from the initial state I^Fq). Once the stirring has 
been stopped, the dominant Hamiltonian in the laboratory frame cor- 
responds to the single-particle motion in the trapping potential 

^trap = Y,( m + P) d m,P dm >P " 
m,p 



With respect to ^ rap , the state ) has the same energy as l^o) 
and the coupling between these two states is thus the dominant es- 
cape route from l^o)- To estimate the corresponding rate, we con- 
centrate on the last part of the adiabatic evolution and we take l^o) 
equal to the Laughlin state, where L = N(N — 1). The coupling ma- 
trix element is then 




Since and l^o) have the same trapping energy, the detuning 

|£!_2 L — Eq\ originates solely from the difference in interaction en- 
ergy. More precisely a lower bound for this detuning is the Laughlin 
gap A Lau ~ 0.1c2- Hence, for N = 10 and C2 = 0.5, the non-resonant 
condition T 1 ^ <C \E}^ — Eq\ is satisfied if the defect amplitude u 
is much less than 10 -2 . For u = 10 -2 , we expect the population in 
|*Fo) t0 decay in a time on the order of 1 /T 1 ^ ~ 20. 

The off-centering defect J%? v expanded at first order in v connects 
|*Fo) on ly with states |¥ a ) whose energy detuning is equal to 1. 
This large detuning is favourable to minimise the departure rate from 
l^o)- Moreover, the influence of this defect fades away together with 
the plug during time evolution. Hence, repeating a similar analysis 
as for Jif u , we eventually find that the condition T a <C \E a — Eq\ is 
safely fulfilled when v < 1 (in units of the trap length), which is not 
a very stringent condition in practice. 
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